Separability criteria for genuine multiparticle entanglement 
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We present a method to derive separability criteria for the different classes of multiparticle en- 
tanglement, especially genuine multiparticle entanglement. The resulting criteria are necessary and 
sufficient for certain families of states. This, for example, completely solves the problem of classifying 
iV-qubit Greenberger-Horne-Zeilinger states mixed with white noise according to their separability 
and entanglement properties. Further, the criteria are superior to all known entanglement criteria 
for many other families; also they allow the detection of bound entanglement. We next demonstrate 
that they are easily implementable in experiments and discuss applications to the decoherence of 
multiparticle entangled states. 
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I. INTRODUCTION 



Entanglement is relevant for many effects in quantum 
optics or condensed matter physics and its characteri- 
zation is of eminent importance for studies in quantum 
information processing [l], Q ■ Concerning entanglement 
between two particles, many questions are still open, but 
there exist at least various criteria which can be used 
to test whether a given quantum state is entangled or 
separable. For more than two particles, however, the 
situation is significantly more complicated, as several in- 
equivalent classes of multiparticle entanglement exist and 
it is difficult to decide to which class a given state be- 
longs. Entanglement witnesses and Bell inequalities can 
sometimes distinguish between the different classes @, Q . 
However, it would be desirable to have useful criteria 
which allow to detect the different classes of multipar- 
tite entanglement directly from a given density matrix; a 
general method to derive such criteria is missing [J]. 

In this paper we present such a systematic way to de- 
velop multiparticle entanglement criteria. The resulting 
criteria solve the separability problem for certain families 
of states (notably, the well-studied A^-qubit GHZ states 
mixed with white noise) and improve known results in 
many other cases. Also, they allow to detect bound en- 
tangled states which are separable under each partition, 
but not fully separable. Moreover, our criteria can eas- 
ily be used in todays experiments and they improve the 
understanding of decoherence in multiparticle quantum 
systems. 

Let us recall the main definitions for multipartite en- 
tanglement. For three particles, a pure state is fully sep- 
arable if it is of the form \ip ls ) = |a)|6)|c) and a mixed 
state is fully separable if it can be written as a convex 



combination of fully separable pure states 

Q ls =J2P^k)(^kl (1) 
k 

where the pk form a probability distribution. A pure 
state is called biseparable if it is separable under some 
bipartition. An example is |"0 bs ) = I a )l0 6c ) where \4> bc ) 
is a possibly entangled state on the particles B and 
C. This state is biseparable under the A|i?C-partition, 
other bipartitions are the B\AC- or C|A_B-partition. A 
mixed state is biseparable if it can be written as g° s = 
SfePfclV'fc S )(V'fc S | where the \ip^ s } might be biseparable un- 
der different partitions. Finally a state is genuine mul- 
tipartite entangled, if it is not biseparable. This class 
of entanglement one usually aims to generate and verify 
in experiments Q and we mainly consider entanglement 
criteria for this type of entanglement. Note that gener- 
alizations and further classifications can be found e.g. in 
Refs. |1 0-13 . 

II. THREE QUBITS 

We explain our main ideas using three qubits, the gen- 
eralization to more particles (or higher dimensions) is 
straightforward and will be discussed later. For a three- 
qubit density matrix g we denote its entries by Qij, where 
1 < i\3 < 8, here and in the following we always use the 
standard product basis {|000), |001), | 111)}. Then we 
have: 

Observation 1. Let g be a biseparable three-qubit 
state. Then its matrix entries fulfill 

1 01,8 1 < V02,2£>7,7 + V £3,306,6 + V 04,405,5 ( 2 ) 

and violation implies genuine three-qubit entanglement. 
Proof. First, note that for two positive linear functions 



2 



f(x) and g(x) the function h = \fjg is concave, that 
is, h[rxi + (1 — r)x2\ > rh(xi) + (1 — r)h(x 2 ) for any 
mixing ratio r [§] . Consequently, the function > /£>2,2£ , 7,7+ 
-y/ 03, 306, 6 + ^04,405,5 ~ 1 01,8 1 is concave in the state, since 
it is a sum of concave functions of the matrix entries 
(the absolute value is convex). So it suffices to prove its 
positivity for pure biseparable states, then mixtures of 
these will inherit the bound. Let \ip) — (cto|0) + ® 

(&oo|00)+& O i|01)+6io|10)+6ii|ll)) be a pure state, which 
is biseparable under the A\BC partition. For that, one 
can directly see that |5i,s| = ^qIaQ^E- For the other 
two bipartitions one finds |5i,s| = -y/03, 3 06,6 and \gi.s\ — 
J 02,207,7) hence, Eq. ([2]) is valid for any pure biseparable 
state, which proves the claim. □ 

This criterion has also been derived in the context of 
quadratic Bell inequalities @, however, our proof is con- 
siderably shorter and, most importantly, it can be gen- 
eralized to derive other characterizations of the different 
entanglement classes. Note that Eq. © is independent 
of the normalization of the state, simplifying many calcu- 
lations below. Eq. © is maximally violated by the GHZ 
state, \GHZ 3 ) = (|000> + |lll>)/%/2. For other states, 
one may first change the local basis (leading, e.g., to the 
criterion \g 2 , 7 \ < ^/5i,i5 8 ,8 + V 03,306, 6 + V 04,405, 5), but 
these will not be considered as independent criteria. 

To discuss the strength of Observation 1, we con- 
sider states which are diagonal in the GHZ basis. This 
basis consists of the eight states \ipi) — (\x1X2X3) ± 
\x1X2X3}) / where Xj,Xj G {0, 1} and Xj 7^ %. States 
which are diagonal in this basis are of the form 
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with real Ai and /Ui, fulfilling Ai = Ag_i for i 
and AT denotes a normalization. We can state: 

Observation 2. For GHZ-diagonal states, the crite- 
rion from Observation 1 constitutes a necessary and suf- 
ficient criterion for genuine multipartite entanglement. 
Proof. The proof is given in the Appendix. □ 

This shows that the criterion of Observation 1 is a 
strong criterion in the vicinity of GHZ states, indeed its 
later generalization solves the problem of classifying N- 
qubit GHZ states mixed with white noise (see Fig. 1). 

It remains to investigate what happens for other states, 
such as the W state, \W 3 ) = (|001) + |010) + |100>)/a/3. 
First, one can apply local unitary operations, before test- 
ing Eq. $2$. This indeed works for the pure W state, but 
one can also derive stronger criteria: 

Observation 3. Any biseparable three-qubit state 




FIG. 1: (color online) The entanglement properties of N- 



qubit GHZ states mixed with white noise, g 



(ghzN) 



(1 



p)\GHZ N )(GHZ N \ + pl/2 N . It was known before that 
these are fully separable iff 1/[1 + 2 (1 ~ JV) ] < p < 1, while 
for smaller p they are inseparable under any partition. Our 
results show, that iff < p < 1/[2(1 — 2~ N )] the states are 
genuinely multipartite entangled. Consequently, in the region 
in between the two bounds the states £A ghzN ' are biseparable 
yet inseparable under any fixed bipartition. 



fulfills 



\Q2,s\ + 1 52,5 1 + 1 53,5 1 < x/01404,4 + ^01, 106,6 + ^/WJoTJi 



+ ^(02,2 



03,3 + 05, 5j 



(4) 



Proof. Again, it suffices to consider pure states. Then, for 
a state which is ^BC-biseparable one sees that | £?2,5 1 = 
-y/0i,i06,6 and 1 03,5 1 = ^/5i,i07,7- Furthermore, one has 
1 02,3 1 < (02,2 + 03,3)/2 which follows already from the 
positivity of the density matrix. Therefore, Eq. (QJ holds 
for the A\BC partition, and similarly one can prove it 
holds for the other two bipartitions. □ 
This observation deserves two comments. First, this 
criterion is quite strong: It detects W states mixed with 
white noise, i.e., 5 (w3) (p) = (1 ~ p)\W 3 ){W 3 \ + pi/8, 
for p < 8/17 sa 0.471 as genuine tripartite entangled, 
while the best known entanglement witness detects it 
only for p < 8/19 « 0.421 [l(J. Second, it should 
be noted that Observation 3 is independent of Obser- 
vation 1: The states p (w3) (p) for p G (0.413; 0.471) are 
directly detected by Observation 3. They are, however, 
not detected by Eq. @, even if local filter operations 

g H> q = T\ ® 7^2 ® -7 r 30-7 r i <8> ^2 ® -?3 are applied with 
arbitrary matrices as can be checked numerically. 

So far, we have only considered criteria for biseparable 
states. Our approach also allows to derive criteria for 
other entanglement classes: 

Observation 4. (i) For fully separable three-qubit 
states, the following inequalities hold: 

i 

\Ql,s\ < (02,2 ' 03,3 ' 04,4 ' 05,5 ' 06,6 ' 07,7) 6 (5) 
|01,s| < (01,1 ' q\a ■ 05,5 ' 06,6 ' 07,7) 6 (6) 

(ii) Eq. ([S]) and Eq. © are connected via the substitu- 
tion 02,203.3 - > 0i.i04,4- Similarly, one obtains new sep- 
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arability criteria from Eq. ([5]) by making the substitu- 
tions 06,607,7 -> 05,508,8, 02,205,5 ~> 01,106,6, 04,407,7 ~> 
03,308,8, 03,305,5 ~» 01,107,7, and 04,406,6 ~» 02,208,8- 

Combining such substitutions, one also obtains new sep- 
arability criteria, e.g. \gi t g\ < (02,2 ■ 03,3 • 05,5 • 08,s) 1/4 - 

(iii) A condition for full separability which is violated in 
the vicinity of a W state is 

|02.3| + |02.5| + 1 03,5 1 < \/01, 104.4 + y/QljQ6fi + ^1,107,7- 

(7) 

(iv) Eq. ((3J) is a necessary and sufficient criterion for full 
separability for GHZ states mixed with white noise. 
Proof. The proof is essentially the same as before, us- 
ing the concavity of more generalized functions [9|. The 
inequalities (|5I7I) are equalities for pure fully separable 
states. The substitutions as in Eq. §Q can be made, 
since 02,203,3 = 01,104,4, etc. holds for any pure fully sep- 
arable state. Concerning (iv), note that Eq. ([5]) detects 
noisy GHZ states for p < 4/5, and this value is known to 
mark the border of the fully separable states Q. □ 

Surprisingly, substitutions as in Eq. do indeed im- 
prove the criterion in some cases. For example, consider 
the family of bound entangled states of Ref. These 
are states as in Eq. ((3]) with Ai = As = /Ji = 1 and 
A 2 = I/A7, A3 = l/Ae, A 4 = I/A5, and fi 2 — £*3 = M4 = 0. 
For A2 • A3 7^ A4 these states are separable under each bi- 
partition, but not fully separable. Their entanglement is 
detected by Eq. ^ or other substitutions. Moreover, as 
one can directly check, for the special case A2 = A3 = A 5 
the inequality in (ii) tolerates significantly more noise 
than the best known witness [flj and gives more signifi- 
cant results for recent experiments. [12J. 

III. MANY QUBITS 

Let us start with introducing a compact notation. 
First, we label the diagonal elements of g by the cor- 
responding product vector in the standard basis. That 
is, if I = (ii,«2) —jijv) is a tuple consisting of N indices 
ik € {0,1} then gi = Qn X! i 2j ... t i N ) is the diagonal entry 
corresponding to \ix,i2, *at)(zi, *2, z jv | - For example, 
for three qubits 0(000) = 0i.i and 0(001) = 02.2 etc. For 
a given / one can define / as the tuple arising from / if 
zeroes and ones are exchanged, e.g., (001) = (110) Fur- 
thermore, let |/| denote the number of ik = 1 in /, then 
£m=« denotes a sum over all / with \I\ = n. 

Second, let a = be a target state and be a 

different state. We abbreviate with O^(g) the sum of 
the absolute values of the off-diagonal elements of g in the 
upper triangle, which correspond to matrix entries where 
a does not vanish. For instance, for the three-qubit GHZ 
state we have £)\ GHZ z) (g) = 1 01,8 1 and Eq. ^ can now be 
conveniently rewritten as 0^ w ^ (g) < X)m=2 y/ 6 (000)61 + 
2"£|I|=1 

The idea behind this notation is to estimate all off- 
diagonal elements similarly as in Observation 1. Ex- 
plicitely, we have for four qubits: 



Observation 5. (i) From the four-qubit GHZ state, 
\GHZ 4 ) = (|0000) + |H11))/V2, a necessary condition 
for biseparability of a general state g is 

^ GHZ ^{Q)<\ £ V^f- (8) 

m={i,2,3> 

This condition is necessary and sufficient for biseparabil- 
ity of GHZ-diagonal states in the sense of Observation 2. 

(ii) From the four-qubit W state, |W 4 ) = (|0001) + 
|0010) + |0100) + |1000))/2, a criterion is derived as 

£)IW>( e ) < J- ^(0000)07 + £ 0/- (9) 
|/|=2 |/| = 1 

(iii) From the four-qubit Dicke state, \D&) = ( 1 001 1) + 
|0101) + 1 1001) + |0110) + 1 1010) + |1100))/V§, a criterion 
is derived as 

O l ° 4) (g) < V£(oooo)0(im)+ £ £ \ZQl0J+2 £ 0I ' 

|J|=1|J|=3 |/|=2 

(10) 

Proof, (i) is proved as in Observations 1 and 2. The 
factor 1/2 takes into account that each possible term 
occurs twice in the sum. (ii) and (iii) follow as in 
Observation 3. Here, estimating an off-diagonal ele- 
ment can be simplified by the following rule: If the 
off-diagonal element rj corresponds to | z 1 «2 * 3 *4 ) ( J 1 J2 J 3 1 
and the state is separable under the A|i3C-D-bipartition, 
one has 77 < y/WJ^7jIJi)BU^hU) wnile one nas V - 
\/Q{i^hu)Q{3ihi 3 iA) for the -4£|CL>-bipartition, etc. 
Further, one needs that for a positive n x n matrix P 
the bound £ 2<j \P tJ \ < ^Tr(P) holds 0. □ 

Again, these criteria improve known conditions: For 
the four-qubit W state mixed with white noise, Eq. © 
detects genuine multipartite entanglement for p < 4/9 ss 
0.444, while the fidelity based witness detects it only for 
p < 4/15 ps 0.267 and the improved witness [Io[ for p < 
16/45 ps 0.356. A four-qubit Dicke state mixed with white 
noise is detected by Eq. ([10]) for p < 8/21 ps 0.381, while 
the best known witness detects it for p > 16/45 ps 0.356 

For arbitrary states similar entanglement criteria can 
be derived as follows: In a given basis and for a fixed par- 
tition, any off-diagonal element can be estimated as in the 
proof of Observation 5. Then, all these estimates can be 
summarized to an estimate of the sum of all off-diagonal 
elements. This might be further improved by considering 
a weighted sum. For instance, for iV-qubit GHZ states, 
the criterion reads D^ GHZn ^ (g) < | S|/|=i yjQi&ji and 
is again necessary and sufficient for GHZ diagonal states 
as the proof of Observation 2 can directly be general- 
ized (see Fig. 1). Further criteria for cluster states or the 
four-qubit singlet state will be presented elsewhere. 
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IV. EXPERIMENTAL CONSEQUENCES 

Obviously, these criteria can be applied to experi- 
ments where the full density matrix has been deter- 
mined [l5[. However, often this can not be done. Still, 
our results may be directly applied. For example, let 
us consider Eq. Q for the detection of entanglement 
around the three-qubit W state. Using the fidelity F — 
Tr(g\W 3 )(W 3 \) one may rewrite Eq. flj) as 

2, s 

F < -(y 'Ql,lQ4,4 + y/Ql,lQ6,6 + \/Ql,lQ7, 7 + £2,2 + 03,3 + 05,5,) • 

The fidelity of the W state can be measured experimen- 
tally with five local measurements [l6j and the diagonal 
elements can also be determined from measurement of 
a z ® it z ® a z , which is already included in the measure- 
ments needed for the fidelity. This shows that Eq. (jJJ) 
(and similarly all other criteria presented) is experimen- 
tally easily testable. For the usual error models in photon 
experiments one can also check that criterion (QJ detects 
entanglement with a higher statistical significance than 
the witness, unless the fidelity is close to one and the 
significance of both methods is high. 



V. DECOHERENCE 

Finally, our results also shed light on the decoher- 
ence of multipartite entanglement. Consider an N-qubit 
GHZ state, influenced by relaxation - the noise that 
is dominant in ion traps On a single qubit, this 

changes the density matrix according to |0)(0| — > |0)(0|, 
|1)(1| -> x|l)(l| + (1 - z)|0)<0|, and (|0)(1| + h.c.) -> 
x 1 ' 2 (|0)(l| + h.c.) (with x — e~ 7 *) and corresponds to 
a coupling to a bath with zero temperature. The to- 
tal density matrix can directly be computed [l8j . result- 
ing in Q12N = x N l 2 for the off-diagonal element and 
Qi = [8\i\ l o+x> I >(l — x) N ~\ I >]/2toT the diagonal elements. 
Here, we have used the same notation as in Observation 
5. 

This state is not diagonal in the GHZ basis, but ap- 
plying on each qubit a filter g t— > J-gJ- with J- = 
a\0){0\ + (l/a)|l)(l| and a 4 = x/(l - x) maps it to a 
state that differs from a GHZ diagonal state only in the 
element Q\ t \. This filtering keeps all entanglement proper- 
ties, but finally Observations 1 and 2 can be used. From 
this one can conclude that GHZ states coupled to a bath 
with zero temperature are genuine multipartite entan- 
gled, if and only if t < - ln[l - (2 N - 1 - l)- 2 ' N ]h. 

VI. CONCLUSION 

We presented a method to derive separability criteria 
for different classes of multipartite entanglement directly 
in terms of density matrix elements. The resulting crite- 
ria are strong and can be used in experiments, as well as 



for the investigation of decoherence. It would be inter- 
esting to use our approach to discriminate between more 
special entanglement classes (such as the W and GHZ 
class for three qubits Q) and to connect it to the quan- 
tification of entanglement with entanglement measures. 

We thank J. Uffink for fruitful discussions. This work 
has been supported by the FWF (START prize) and the 
EU (OLAQUI, QICS, SCALA). MPS acknowledges the 
hospitality of the Centre for Time, University of Sydney. 

Note added: Half a year after submission of our 
manuscript to the arxiv, a preprint [l^| has appeared, in 
which criteria for multipartite entanglement have been 
presented. This method is analog to ours, by u sing esti- 
mates for off-diagonal terms [cf. Eq. (I) in Ref. [19( with 
the proof of our Observation 5] and using convexity argu- 
ments. Consequently, the obtained separability criterion 
used for N-qubit GHZ states [Eq. (II) in Ref. 119]] is the 
same as our criterion in Eqs. ([2j [8]) combined with lo- 
cal filtering operations, and the criterion for W states 
[Eq. (Ill) in Ref. [19]] is for three qubits the same as 
our Eq. (jj), while for four qubits it is weaker than our 
Eq. ©. 

Appendix 

Here, we prove Observation 2. Since g>( dla ) > Owe have 
\m\ < Aj and we can assume that Ai > A2 > A3 > A4 as 
one can achieve that by a local change of the basis. Then, 
Eq. ([2]) reads \m\ < A2 + A3 + A4, and we will show that 
if this holds, a decomposition into biseparable states can 
be found. Note that due to the ordering of the Ai other 
conditions for biseparability (e.g. \^\ < Ai + A3 + A4) 
can then never be violated. 

Let us define the unnormalized state g^ 12 \X) with Ai = 
fix = A2 = fi2 = A, while all other matrix entries vanish. 
This state is A£?|C-biseparable, since it can be written 
as 

gW(\) =2\(\ X + )(x + \ab®\v + )(v + \c+ 

+ lx~Xx~U B ®brXTk) (11) 

with | X ± ) = (|00) ± |11»/V2 and 1^) = (|0) ± |1»/V2. 
Analogously, one can consider states g^ for any k, I = 
1, 4 with k 7^ I and finds that they are also biseparable, 
as one only has to permute or flip some qubits. 

(i) First, we consider the extremal case when ^ = A; 
for all i and by assumption the separability condition 
implies that we have [ii — \ < X^fe^i Afe, where for the 
index 1 < k < 4. If Ai = A2 + A3 + A4 we can directly 
write £i( dla ) = J2k=2 3 4 6 l< ' lfe '(Afc) hence p( dla ) is bisepara- 
ble. Otherwise, the idea is to write 

P (dia) = E e {lk) (x k ) + e M (12) 

fc=2,3,4 

for some parameters Xfc such that the rest g^ (which is 

(r) 

then characterized by parameters X k ) fulfills two con- 
ditions. Its first and last column and row should vanish 
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(f?ii = ^i = 0) and it should still fulfill all biseparabil- 

ity conditions (e.g. < A3 + A^). Then, g^ can be 
iteratively further decomposed and finally a decomposi- 
tion of £i( dla ) into biseparable states can be found. 

The idea is to choose the \ k ,k = 2,3,4 as equal as 

possible (they have to fulfill Aj, < A&), but monotoni- 
cally decreasing. For that, we define 04 := A2 + A3 + 

(r) (r) (r) 

A4 — Ai = A 2 + A 3 + A 4 > and then recursively 

A 4 = min{A4, 04/3}, then 03 = 0:4 — A 4 and then 

A3 = min{A 3 ,ck3/2} and finally ev 2 — a 3 — A3 and 

A 2 r) = min{A 2 ,a 2 }. Then g^ = Ek=2, 3 A 9 {lk) (A* - 

Ai r) ) + g (l) with A 2 r) > A^ r) > \%\ Then we cannot have 

that both A4 = A4 and A3 = A3, because if this were 
true, then from the definition of 0:4 it would follow that 

A 2 r) = A 2 - Ai < 0. So we have A 3 r) = ev 3 /2 (A 4 r) = a 4 /3 
also implies A3 = 03 /2), which due to the ordering of 
the Ai implies A 2 r) = a 2 = A 3 r) . So A 2 r) < + A 4 r) and 
one can decompose g^ further into g^ 23 ^ and g( 24 ^ and a 

(rr) (rr) 

remaining term with \\ — A 2 — etc. Of course, for 
the case of three qubits one may also write down suitable 
values for the A^ directly, but the previous scheme can 
directly be extended to more qubits. 

(ii) Secondly, for < fa < Ai, and where again 
Ai < A 2 + A 3 + A 4 , we first consider the states g( kl K 
Their non-zero matrix elements obey fa — Ai, but, ap- 
plying with some probability locally conjugate random 



phases (e.g, |1) 2 i-> e**|l) 2 and |1) 3 n- e _i *|l) 3 ) to these 
states decreases the values of the fa (in this example for 
i = 2,3). Therefore, one can for a given g( k[S> decrease 
the values of fa arbitrarily by local operations (in the ex- 
ample we can decrease e.g. the value of fa for g/ 12 * 1 or fa 
for g^ 34 ^), and the resulting states must be biseparable. 
Consequently, a given p( dla ) with Ai < A 2 + A 3 + A4 can 
be decomposed into biseparable states as in (i). 

(iii) Further, it may happen that for a given p( dla ) one 
has < fa < A 2 + A 3 + A 4 but Ai > A 2 + A 3 + A 4 . Then 
we consider g which is obtained from p( dla ) by setting 
Ai = max{/ii, A 2 }. Then, g is biseparable according to 
(ii), and £A dla ' is obtained from g by mixing with the 
fully separable state |000)(000| + |111)(111|, hence it is 
biseparable. 

(iv) The previous arguments prove the claim if all 
fa > 0. If some fa are negative, one can prove it as fol- 
lows: Let g be a GHZ diagonal state, with some fa < 0, 
which fulfills the condition of biseparability. The state g 
which arises from g when all fa are replaced by \fa\ ful- 
fills the same condition, and is biseparable due to points 
(i)-(iii). It can be decomposed into several g^ kl \ maybe 
in some of them we have ) < Xi{g^) according 
to points (ii) and (iii). Nevertheless, we can built out 
of this decomposition of g a decomposition of g, if we 
flip the signs of all the fa(g^) appropriately. An arbi- 
trary flipping of the signs of the fa of a given g( kl > can 
be done for each k, I by local operations, hence g is also 
biseparable. □ 
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